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We present a modification of Matrix Product State time evolution to simulate the propagation of 
signal fronts on infinite one-dimensional systems. We restrict the calculation to a window moving 
along with a signal, which by the Lieb-Robinson bound is contained within a light cone. Signal 
fronts can be studied unperturbed for much longer times than on finite systems. Entanglement 
inside the window is naturally small, greatly lowering computational effort. We calculate the time 
evolution of the transverse field Ising (TFI) model and of the XXZ spin-half antiferromagnet after 
local quantum quenches. In both models, we observe distinct magnetization plateaus. In the TFI 
model, their dynamical scaling follows power laws to high precision. 



Propagation of signals in one-dimensional (ID) 
strongly interacting systems has been of longstanding 
general interest in both condensed matter and quantum- 
computational physics. The phenomenon provides a ba- 
sis for coherent information transfer via quantum wires. 
A signal can be created, e.g., as a local excitation from 
a stationary state, or as a domain wall or a topologi- 
cal excitation [U[2]. In recent years, the nonequilibrium 
time evolution of such signals after global jSHTT] and lo- 
cal p!2H2Q] quantum quenches has been the subject of 
intense theoretical interest. Often hard to pursue by an- 
alytical methods, many studies have become feasible in 
ID by numerical methods [2TJ[22 based on Matrix Prod- 
uct State (MPS) [23l [24], and also by tensor network 
approaches [25 . 

If the initial state is prepared within a finite system, 
boundary effects such as Friedel oscillations interfere with 
a passing signal. System boundaries also limit the time 
span for tracing the signal, before it is reflected at the 
boundaries in a non-trivial manner. The maximum time 
is even more severely restricted by entanglement which 
develops across the system and which requires a compu- 
tational effort that can drastically increase with time [26] . 
This has greatly hampered the analysis of large time 
asymptotic behavior [8l |9l |27l [SHj . Boundary effects do 
not appear in infinite systems, for which the ground state 
and its time evolution can be calculated with MPS meth- 
ods [28H3T] . However, these methods require complete 
translational invariance and can therefore not be applied 
to study signal propagation. 

In this letter we present a method to simulate the prop- 
agation of local signals on an infinite chain using MPS 
time evolution, without any finite size effects distorting 
the signal front. We study the time evolution of the 
Transverse Field Ising (TFI) model [3H3 HH HI |32H3Il 
and of the XXZ spin \ chain (SHIIl [HHIIl [38] after a 
local quantum quench up to large times. We find dis- 
tinct magnetization plateaus [HI [38] , reminiscent of the 
plateaus observed in equilibrium frustrated spin chains 



as a function of magnetic field j 39ti46] . In our case the 
plateaus develop over time in the longitudinal magnetiza- 
tion close to the signal front. Their heights and positions 
scale algebraically with time in the TFI model to high 
precision. 

For our method we consider a general spin chain of infi- 
nite size, with finite range interactions, initially prepared 
in a state, such as the ground state, which is translation- 
ally invariant for sites x > x^io the right of some site xq. 
At time zero, the system is excited by a quantum quench 
like one or more spin flips at sites x < or a modifi- 
cation of the Hamiltonian at x < xq. For finite range 
interactions it is known from the Lieb-Robinson bound 
[47] that wave fronts generated by such quenches can 
only propagate with a characteristic maximum velocity 
^max, i-e. within a "light cone" even in a non-relativistic 
system, as recently also seen experimentally [Ij. Any 
correlations beyond the light cone are exponentially sup- 
pressed at increasing distance with a length scale of at 
most the correlation length [47 . 

We divide the system into three parts, namely the co- 
moving window (CMW) which moves towards the right 
with the wave front, and two half- infinite regions, a uni- 
form one in front (i.e. to the right) of the window, and 
an arbitrary one to the rear. The window is chosen wide 
enough to contain the complete signal front, including 
the exponentially damped part to the right of the main 
front, to high precision. The signal will therefore not 
affect the system to the right of the window. Likewise, 
when the window moves with Vmaxi modifications in the 
rear part will not affect the CMW and need not be cal- 
culated. Since bipartite entanglement [48H5Q] spreads at 
most with Vmaxi the bipartite entanglement entropy is 
significantly lower around the wave front than in the 
bulk, allowing for low computational effort when using 
the CMW. 

We label the sites according to the division we have in- 
troduced: left side j < CMW £ + 1 < j < r, and right 
side j > r + 1. We consider a Hamiltonian with near- 
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Figure 1: (color online) Time evolution of magnetizations 
{S^{t)) and {S^{t)) of the TFI in a Comoving Window af- 
ter a Jordan- Wigner excitation, at h — 0.45. Inset: Time 
evolution after a single spin flip of without window move- 
ment. Color scale spans the data, from 0.074 (blue) to 0.4063 
(red). 

est neighbor couplings H = hjj-^i^ which subdivides 
correspondingly into 

H = H^^ ^£,£+1 + + K,r+1 + • (1) 

The wave function in the basis becomes 

^(. . . Sl-iSl I 5£+lS£+2 ...Sr\ Sr+lSr+2 • • • ; • (2) 

where Sj labels the spin at site j. It can be approximated 
very closely in the form of an infinite matrix product [23 

"(3) 

where L^3<^ are left-orthogonal matrices in the left part, 
j^si+i<j<r- ^s£+i<j<r right-orthogonal 
matrices inside the CMW, respectively, R^^+^<j are right- 
orthogonal matrices in the right part and A^-^-^ are di- 
agonal matrices containing the Schmidt values of a bi- 
partition at bond {k^k + 1). For a finite system, the 
left and right ends of (|3| are terminated by contrac- 
tions with boundary vectors; we however consider the in- 
finite size limit. In standard TEBD notation [22 , L^^ = 
A^-ir^^A^^- = X^-^r'^,B'^ = r^^A^ and R'^ = T'^X^. 
The matrices R^^ are obtained as the stationary state of 
iTEBD[28], thus there are only two independent matri- 
ces, and R^^+^ R^^ . 

Let us consider one step of the time evolution. Us- 
ing Suzuki- Trotter decomposition [51 , we apply adaptive 
time evolution with respect to Hu + /^r,r+i + inside 
and to the right of the CMW, and, e.g., non-adaptive 
time evolution with respect to Hi, + hi^-^i to the left of 
the CMW. 

Adaptive time evolution inside the CMW can be per- 
formed either by tDMRG [2l] or TEBD [22], using local 
operators Uj j_^i{r) = e~'^^^o,3+i with finite time steps 
r. To the right side of the CMW, the evolution of 
ig performed by iTEBD [28ll52]. 
Special attention has to be paid to the operation of 
(r) at the junction between the CMW and the right 
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Figure 2: (color online) Magnetization {Sx{t)) at times 190 < 
t < 500, showing magnetization plateaus. Inset: {S^{t)) (red 
crosses) vs. position inside the window at t = 100 (lower line) 
and t = 300 together with the exact solution (blue circles, 
coinciding with data). 

side. We explain our procedure in terms of tDMRG rep- 
resentations. After a sweeping process |21j to the right 
boundary of the CMW {j = r), and the first iTEBD 
half-step (on f^^^ ^i^q following ma- 

trix product 

. . . A''+' A''+^ ...A'^ X'^R'^+' R'^+^ .... (4) 

We then operate with 'U^^^+i(r) on t/j^r-Sr^+i ^ A^'^X^R^'^+^ 
and apply singular value decomposition (SVD including 
truncation) on the result, which yields 

with cj^ containing the updated Schmidt values. 

After the second iTEBD half-step (on r^-+2 A^+^T^^+i), 
the matrix T/^'^+i will differ from the new matrix R^^+'^ 
by some unitary MPS gauge transformation x, such that 
ysr^i _ ^ i^«r-+i ^ since there is negligible influence of 
the signal around the right boundary by construction 
(see also results below). Using right orthogonality [23] 
R'R'^ = 1, X is obtained as x = V'-+^ R'-^^K 

We arrive at 

^5£+i^5.+2 A'^-^V'-uj^x ^ 

where the right side is again purely represented by ma- 
trices R^K Regular tDMRG can now be continued by 
forming '^^^-i^^ = A^''-^U^''ijf x and sweeping back to 
the left of the CMW, which completes one time step of 
the adaptive update. Note that this protocol can also be 
performed using TEBD. 

Non-adaptive time evolution to the left of the CMW 
can be performed similar to ref. 53, Consider ini- 
tially in a renormalized form with respect to the 
MPS matrices L^^ , obtained, e.g., from an initial DMRG 
calculation [29l [SQ. The matrix is thus defined 

in a basis of the auxiliary space connecting L^^ 

and A^^+i . We include in DMRG fashion to form 

= ^L^ + ^£/+i in the block spin basis {\a£S£-^i)}. 
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Figure 3: (color online) Scaling of {Sx{t)) vs. (x — xo)/t. In- 
sets: Algebraic scaling fits (blue lines) for the first (1) and 
second (2) magnetization plateau (red crosses). Left inset: 
position, right inset: height (S^). 

j[si+i -g updated by [/^'^+^(r) = exp ( - itHl^i+i) 
to get 

This update is non-adaptive because possible changes of 
the matrices are compressed to A^^+i. The CMW 
however, when moving with Vmaxi will not be affected. 
The non-adaptive update can also be done by means of 
the renormalized imaginary-time transfer matrix, as in 
finite-temperature DMRG [54]. An adaptive update with 
gauge transformation could be used as well, when the 
initial state is translationally invariant to the left. An 
even simpler alternative is to set = 0, disconnecting 

the left part, with similar results. 

The window is moved by Vmax when known before- 
hand, or the movement is triggered by the change of 
some local observable or the bipartite entanglement at 
a sufficient distance from the right boundary. A shift of 
the CMW is done by including R^""^^ as a new rightmost 
matrix into the CMW, and by absorbing the leftmost 
matrix A^^+i into the left side by using it to renormalize 
^L^+i- Similar operations apply for an alternative shift 
by a whole iTEBD unit cell. 

Results. The spin-^ Transverse Field Ising model 
[2 El [3 IpM] on an^ infinite chain defined by 
H = — Yljp^^^+i ~ '^j can be diagonalized ex- 
actly [551 [56] by means of a Jordan- Wigner (JW) trans- 
formation [57] , but the time evolution of observables is in 
general very expensive to calculate [33 . At h < he = 0.5 
there is long range order in . We prepared the system 
in the symmetry broken ground state |0) with (Sj) > 
using iDMRG [29l |30]. As a local quantum quench 
we employed an initial JW excitation (cj.^ -\- Cx^) \0) = 

U^<x, (-2^z") (^^0 + ^xo) |0) at site xo inside the win- 
dow, where c^,c are JW fermion creation and annihila- 
tion operators [57] • It corresponds to a spin flip in z- 
direction at site xq and a domain wall in x-direction be- 



tween sites Xq — 1 and xq. We present data at = 0.45, 
using r = 0.002 and maximum matrix size 26. 

The time evolution inside the CMW (Fig. [T]) shows 
that boundary effects are indeed removed at both ends 
of the CMW. Window movement was triggered by en- 
tanglement, resulting in a window velocity of 0.455(4), 
equal to the maxmimum group velocity Vmax = h of the 
TFI [3 . When the window is not moved (Fig. [l] in- 
set), an impinging signal is absorbed by both boundaries 
temporarily, though reflections occur eventually. The en- 
tanglement entropy (not displayed) shows reflections at 
about t = 70 at both boundaries. 

The agreement of the transverse magnetization in 
Fig. |2|inset) with exact results which we calculated for 
the JW excitation is excellent, even at very large times, 
with relative deviations of less than O(10~^) beyond the 
rear boundary. The transverse magnetization develops 
two shallow domain walls after the quench, with an over- 
all scaling {S^ — SUt)) ^ 1/t away from the signal front, 
as expected [33l|58.- Close to the front, there are oscil- 
lations. The difference of the ground state magnetiza- 
tion and the first minimum numerically fits a power law 

^ ^-1.127(8) _ ^-9/8^ 

For the longitudinal magnetization (6'^(t)), analytical 
calculations are much more difficult, and results have 
been available in the literature only very recently [3 . We 
find that shoulders appear during the time evolution, 
which develop into distinct magnetization plateaus, with 
a relative variation in {S^{t)) of less than 0.1% over 4 
sites for the rightmost plateau in Fig. [2j The plateaus 
slowly widen over time. Analoguous steps have previ- 
ously been observed in the XX chain [18 , and with (S^) 
in the paramagnetic phase of the TFI at = 5 \3E\. 
Fig. [3] shows that the magnetizations at different times 
can be scaled to a single function {S^{t)) = f{{x — xo)/t) 
away from the front, similar to the behavior of {S^{t)) in 
Ref. |35| The plateaus close to the signal front are excel- 
lently described by a power law {S^{t)) = m-\-at^, where 
m = 0.4026(2), b = -0.3497(6) for the first plateau 
and m = 0.3982(4), b = -0.362(1) for the second one, 
with m = 0.406269 for the ground state. Their posi- 
tions p obey {p — xo)/t = v — ct^ with v = 0.4500(3), 
c = 1.52(2), d = -0.675(4), resp. v = 0.4503(4), 
c = 2.51(3), d = —0.666(4). The exponents d are very 
close to —2/3. We note that the velocities {p — j)/t con- 
verge to V = h = Vmax to high precision making the 
slope in Fig. [sj^left inset) sharper for larger times. The 
line of crossing points visible in the figure at large {S^) 
thus also converges to v = h. However, since d > —1^ 
the positions of the plateaus move away from the signal 
front with a power law ~ t^~^^ ^ t^^^ . The exponents are 
virtually unchanged within error bars for different fields 
h = 0.2,0.3,0.4 < he (see also [33 ) and for a different 
quench without the local flip. 

Finally, we investigated the XXZ antiferromagnet [8]- 

m MM EH], ^ = -J"^E,(^i^i+i + - 
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Figure 4: (color online) Time evolution of the bipartite entan- 
glement entropy and of {{ — lySj) in the XXZ antiferromagnet 
after a local spin flip at J^^ = 1 and — —3. 
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Figure 5: (color online) Magnetization (( — 1)^5'|) in the XXZ 
model at times 20 < t < 200, showing developing plateaus. 
Inset: data plotted vs. position inside the CMW, for times 
t > 60, with the last 2 plateaus. 

Sj Sj_^i, using r = 0.01 and matrix dimension 
100. We prepared the system in the symmetry broken 
infinite size ground state at J^^ = 1,J^ = —3, with 
(— 1)-^ (Sj) > 0. At time we flipped a spin at site xq, 
thereby producing two domain walls (spinous). During 
the time evolution, these domain walls spread in opposite 
directions (Fig. [4]), forming a widening region with neg- 
ative (-1)^' (^l) in between. We find Vmax = 1.945(2) 
from the evolution of (S^). Close to the initial excita- 
tion, the bipartite entanglement starts to grow consider- 
ably but within the CMW it stays small at later times. 
At large times t, distinct magnetization plateaus develop 
again, as shown in Fig. [5j with a relative variation of less 
than 0.1% over 6 sites for the rightmost plateau in the 
figure. Over time, they widen and move away slowly from 
the signal front (Fig. [sj inset). The plateaus do not fol- 
low power laws similar to the TFT The scaling behavior 
in the XXZ model [8, 9 is currently under investigation. 

Conclusions. We have introduced a comoving win- 
dow, employing an easy-to-implement gauge transforma- 
tion in order to connect the CMW to its surroundings. 
We have shown that the CMW can follow the propaga- 
tion of a signal front on an infinite size system unim- 
peded for very long simulation times. We note that even 



when the window is not moved, local signals can be sim- 
ulated on the background of an infinite system, without 
perturbations which would otherwise emanate from the 
boundary. In this case the signal is temporarily absorbed 
by the boundary, though simulation times are eventually 
limited by the size of the window. 

The CMW method is not restricted to evolution of ex- 
citations under uniform Hamiltonians. For example, the 
AKLT model [59 with inhomogeneous bond interactions 
or ID quantum systems under exponential or hyperbolic 
deformation [6Q[|6T] have uniform ground states, whereas 
the Hamiltonian is not uniform. 

The CMW enabled us to calculate the relaxation af- 
ter local quantum quenches in the Transverse Field Ising 
(TFI) model and in the XXZ antiferromagnet with low 
computational cost up to very large times. We have 
found distinct dynamical magnetization plateaus close to 
the signal front in both models. Their scaling is algebraic 
and well described in the TFI, whereas in the XXZ model 
it is more complicated and under current investigation. 

Note added. During the preparation of this letter we 
learned that work on related approaches was being com- 
pleted by two other groups [62l |63] . 
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